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Abstract 

In this paper we analyze a time-reversal experiment in a random underwater acoustic channel. 
In this kind of waveguide with semi-infinite cross section a propagating field can be decomposed 
over three kinds of modes: the propagating modes, the radiating modes and the evanescent modes. 
Using an asymptotic analysis based on a separation of scales technique we derive the asymptotic 
form of the the coupled mode power equation for the propagating modes. This approximation is 
used to compute the transverse profile of the refocused field and show that random inhomogeneities 
inside the waveguide deteriorate the spatial refocusing. This result, in an underwater acoustic 
channel context, is in contradiction with the classical results about time-reversal experiment in 
other configurations, for which randomness in the propagation medium enhances the refocusing. 
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Introduction 

The time-reversal experiments of M. Fink and his group in Paris have attracted considerable attention 
because of the surprising effect of enhanced spatial focusing and time compression in random media. 
The refocusing properties have numerous interesting applications, in detection, destruction of kidney 
stones, and wireless communication for instance. This experiment is in two steps. In the first step 
(see Figure [l] (a)), a source sends a pulse into a medium. The wave propagates and is recorded by 
a device called a time-reversal mirror. A time-reversal mirror is a device that can receive a signal, 
record it, and resend it time-reversed into the medium. In the second step (see Figure[l](fe)), the wave 
emitted by the time-reversal mirror has the property of refocusing near the original source location. 
However, surprisingly, it has been observed that random inhomogeneities enhance refocusing jTJIE]- 

Time-reversal refocusing in one-dimensional propagation media is carried out in \Q llUj . In three- 
dimensional randomly layered media [9] , in the paraxial approximation [31 [2T] , and in random 
waveguides [TCTl IT!?] , it has been shown that the focal spot can be smaller than the Rayleigh resolution 
formula XL/ D (where A is the carrier wavelength, L is the propagation distance, and D is the mirror 
diameter), but the focal spot is still larger than the diffraction limit A/2. Moreover, in [TJ] the author 
propose a setup in a waveguide with random perturbations in the vicinity of the source in order to 
obtain a superresolution effect, that is to refocus beyond the diffraction limit with a far-field time- 
reversal mirror. The setup was inspired from |20j describing this superresolution effect experimentally. 

Time reversal of a broadband pulse, in the case of a waveguide with a bounded cross-section and 
Dirichlet boundary conditions, is carried out in [T^] and Chapter 20]. In underwater acoustics 
the waveguide model (see Figure [T]) has a semi-infinite cross section, and therefore a wave field can 
be decomposed into three kinds of modes: the propagating modes, the evanescent modes, and the 
radiating modes. Consequently, waveguides with bounded cross-sections does not take into account 
radiation losses which can be encountered in underwater acoustics. Moreover, the refocusing enhance- 
ment of the time-reversal experiment in random waveguides is closely related to the propagation of 
the propagating modes power in the random medium. For waveguides with bounded cross-sections 
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Figure 1: Representation of the time-reversal experiment. In (a) we represent the first step of the 
experiment, and in (6) we represent the second step of the experiment. 



the total propagating mode power is conserved and even uniformly distributed after long-range prop- 
agations |10(ll2j. that is why refocusing enhancement can be observed. In underwater acoustics, with 
a far-field time- reversal mirror only the propagating modes contribute to the focal spot, and it has 
been shown [T^l [T7] that the presence of radiating modes produces an effective dissipation on the 
propagating modes, which cannot be fully compensated by time reversal. As a result, in contrast with 
all the results listed above regarding improvement of the time-reversal focusing thanks to the ran- 
dom perturbations in the propagation medium, we show in this paper that random inhomogencitics 
in a context of underwater acoustics deteriorate the refocusing property. The deterioration of the 
refocusing property for the time-reversal experiment has already been observed in [Tj. However, in 
there context the authors considered a thermo-viscous wave model to incorporate viscosity effects in 
wave propagation. In our context, the deterioration of the refocusing is only due to the random in- 
homogeneities and the geometry of the propagation media. Therefore, it is interesting to understand 
these effects and to characterize them on the focal spot resulting from the time-reversal experiment. 
The main result of this paper is the analysis of the radiation losses on the refocused wave in the 
time-reversal experiment in a context of underwater acoustics. In Theorem |3.7| and Theorem |3.8| 
we show that the radiative losses affect the quality of the time-reversal refocusing in two different 
ways. First, as expected according to the results obtained in |15j . the amplitude of the refocused 
wave decays exponentially with the propagation distance. Second, the width of the main focal spot 
increases and converges to an asymptotic value, which is larger than the diffraction limit X oc /{29) 



obtained in Proposition 3.1 (where A oc is the carrier wavelength in the ocean section with index of 
refraction ni, and 9 = \JT— l/ n i)- 

The organization of this paper is as follows. In Section [l] we introduce the underwater waveguide 
model studied in detail in |24j . and in Section [2] we present the mode decomposition associated to this 
model and we derive the coupled mode equations. In Section|3]we study the time-reversal experiment. 
We describe in a simple way the refocused transverse profile in terms of the solution of the continuous 
diffusive model obtained in [15 , and describing the mode-power coupling between the propagating 
and radiating modes. Thanks to this representation we show that the quality of the time-reversal 
refocusing is deteriorated by the radiative losses in the ocean bottom. 

1 Waveguide Model 

We consider a two-dimensional linear acoustic waveguide model. The conservation equations of mass 
and linear momentum arc given by 



,3u „ „ 

p(x,z)~ + Vp = F q , 

1 9p _p 
K(x, z) dt 



(1) 



where p is the acoustic pressure, u is the acoustic velocity, p is the density of the medium, K is the 
bulk modulus, and the source is modeled by the forcing term F q (t, x, z) given by 



F e {t,x,z) =y e (t,x)S{z-L s )e z . 




Figure 2: Illustration of the shallow-water waveguide model. 



The third coordinate z represents the propagation axis along the waveguide. The transverse section 
of the waveguide is the semi- infinite interval [0,+oo), and x € [Q,+oo) represents the transverse 
coordinate. Here, F e represents a point source localized at z — Ls, pointing in the z-direction, 
with temporal and transverse profile given by ^ e (t,x). Let d > be the bottom of the underwater 
waveguide, the medium parameters are given by 

i (n 2 (x) + y/eV(x,z)) if xe[0,d], ze[0,L/e] 

x G [0, +oo), z e (-oo, 0) U (L/e, +oo) 



K(x,z) I j^n 2 {x) if 

p{x, z) — p if x e [0, +oo), z € R, 



or 

x £ (d, +oo), z £ (— oo, +oo). 



and where V is a stochastic process describing the random perturbation of the propagation medium 
(see Figure |2]for an illustration of the random underwater waveguide model). In this paper we consider 
the Pekeris waveguide model. This kind of model has been studied for half a century [22] and in this 
model the index of refraction n(x) is given by 



n{x) 



ni > 1 if x € [0, d) 
1 if x e [d, +oo). 



The Pekeris profile models an ocean with a constant sound speed, and where d represents the ocean 
depth. Conditions corresponding to the Pekeris model can be found during the winter in Earth's mid 
latitudes and in water shallower than about 30 meters. 

The perturbation V is assumed to be a continuous real-valued zero-mean stationary stochastic 
process with 0-mixing properties |19j . More precisely, let 

J~u = Fo,u = c(V r (x, z), x G [0, d], < z < u) and F u ^oo = a(V{x, z), x £ [Q,d], ,u<z), 

we assume that 

sup \P(A\B) -P(A)| < <j>{v). 

This ^mixing property describes the decorrclating behavior of the random perturbation V through 
the nonnegative function (f> € n L 1 ' 2 (K) characterizing the decorrelation speed. Throughout 

this paper, for the sake of simplicity and for explicit computations in Section 3.9.2| and Section 3.9.3 
we assume that the autocorrelation function of the random perturbation V is given by 

E[V(x 1 ,z 1 )V(x 2 ,z 2 )} = l0 (x 1 ,x 2 )e~ a ^- z ^ V{z 1 ,z 2 ,x 1> x 2 ) € [0,+oo) 2 x [0,d] 2 . 

From the conservation equations (JT|), we derive the wave equation for the pressure field, 

A 1 d 2 p _ v pe 

c(x, z) 2 dt 2 
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where c(x,z) — y/ K(x, z) / p(x, z) is the sound speed profile, A = d 2 + d 2 , and c — y K / p. In 
underwater acoustics it is natural to use a pressure-release condition since the density of air is very 
small compared to the density of water. As a result, the pressure is very weak outside the waveguide, 
and by continuity, the pressure at the free surface x = is zero. This consideration leads us to 
consider the Dirichlet boundary conditions 

p(t,0,z) = V(t, z) e [0, +00) x K. 

To study the focusing property of the time-reversal experiment we need to understand the wave 
propagation in the random medium. To do that, we use a separation of scale technique introduced by 
G. Papanicolaou and his coauthors in [2] for instance. The important scale parameters in our problem 
are: the wavelength, the correlation length and the standard deviation of the medium inhomogeneities, 
the propagation distance, and bandwidth of the pulse. This last scale parameter plays a key role in 
the statistical stability of the time-reversal experiment. 

First of all, the model of wave propagation considered in this paper is a linear models, so that the 
pressure p(t, x, z) can be expressed as the superposition of monochromatic waves by taking its Fourier 
transform. Here, the Fourier transform and the inverse Fourier transform, with respect to time, are 
defined by 

/H = J f(t)e^dt, f(t) = i f f(u>)e- iut du. 

As a result, in the half-space z > L$ (resp., z < Ls), we get that p\ui, x, z) satisfies the time-harmonic 
wave equation without source term 

d 2 p(uj,x,z) + d%p(uj,x,z) + k 2 (uj)n 2 (x)p(uj,x,z) + y/eV(x, z)p(u, x, (z)l[o,L/e] 0) = 0, ( 2 ) 

where k(uS) — ui/c is the wavenumber, and with Dirichlet boundary conditions p\u, 0, z) = Vz. The 
source term does not appear in ^ but induces the following jump conditions for the pressure field 
across the plane z — Ls 



(3) 



p{w,x,L^)-p{u),x,L s ) = V(u,x), 
d z p(u},x,L"£) - d z p(u,x,Lg) = 0. 

To study ([2]), we consider this equation as an operational differential equation 

d 2 

■j-jp(uj,.,z) + R(uj)(p(uj, .,z)) + y/eV(-,z)p(uj, ., z)l [0 ^ d] (-)l [0 , L/e] (z) = 

in H = L 2 (0,+oo), where R(uj) is an unbounded operator on H with domain 

V(R(u)) = ff 1 (0,+oo)n J ff 2 (0,+oo), 

and defined by 

R(u)(y) = ^y + k 2 (u,)n 2 (x)y WyeV(R(cj)). (4) 

In the next section we introduce the spectral decomposition of the operator R(oj) [21] . This decom- 
position we will be used in what follows to decompose the field p(ui,x,z) and then understand the 
stochastic effects undergoes during the propagation. 

1.1 Spectral Decomposition in Unperturbed Waveguides 

The spectral analysis of the self-adjoint Pekeris operator Q is carried out in [23]. To use this spectral 
decomposition, we are interested in solutions of Q such that 

p(u>, ., .)l(i s ,+oo) (2) G C° ((L s , +00), -ffo(0, +00) n H 2 (0, +co)) n C 2 ((L S , +00), ff) , 

p{cj, ., .)i ( _ 00 , Ls) (z) e c° ((-00, l s ),h£(o, +00) n h 2 (o, +00)) n c 2 ((-oo, L s ), h) . 

According to [23] , the spectrum of the unbounded operator Q is given by 
Sp(R(u)) - (-oo,fc 2 ( W )] uj^M,...,/? 2 ^)}. 
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The continuous part of the spectrum comes from the fact that our waveguide model is semi-infinite. 
For the discrete part, the modal wavenumber f3j(uj) are positive and 

fc 2 ( W ) < $r (w) M < • • • < < nlk\u). 

Regarding the spectral decomposition, there exists a resolution of the identity II W of R(uS) such that 
My € H,Vr e E, 



IL UJ (r,+oo)(y)(x) = (y, ^-(w, .)) H <j>j(u), x)l (j%+oo) (/3j(w) 2 ) 
rfc 2 (") 

- / (y,</!> 7 (a;,.)) J j07(w,x)d7l(_ oo>fc 2( w ))(r), 



and Vy e V(R{u)), Vr e 



AT(uO 



n«„(r J +oo)(iZ(w)(y))(x) = J] & 

fc 2 (w) 

Let us describe more closely the discrete and the continuous part of the decompositions. 

Discrete part of the decomposition Wj € {l, . . . , ./V(w)}, the jth eigenvector is given in [55] by 

where 



Aj(oj) sin(aj(Lj)x/d) if < x < d 

A J (w)sin(CT J (w))e-^( tJ ) 21 ^ if d < x, 



and 



Aj(u;) 



\ 



2/d 



sin 2 (tTj sin(2crj (cj)) 



1 



O (<•>) 2o-j (w) 

Here, cr 1 (a;), . . . ,<Jiv( w )(w) are the solutions on (0, nik(u>)d9) of the following equation, 

tan(y) - --j=====, 



(•5) 



(6) 



and such that < o\(uS) < ■ ■ ■ < ct n ^(uj) < nik(uj)d9, with 6 = ^/l — l/n\. This last equation ad- 
mits exactly one solution over each interval of the form (n/2+(j — ir/2+jn) for j 6 {1, . . . , iV(w)}, 
where 

7T 

and [•] stands for the integer part. From (|6]), we have the following results [TB] which are used to 
obtain the main result of this paper in Section pHO] 

Lemma 1.1 Let a S (1/3, 1), we have as N(us) — s- +00 

sup |oj+iH - crjM - tt| = O (jVH5"5 a ) . 

je{i,...,Ar(w)-[iV(a»)«]-i} v y 

sup |o-i+a(w) - 2o- J+ i(w) + o-i(w))| = O (N^) 1 ' 30 ) . 

je{i,...,Jv(u)-[jv(w)<»]-2} 
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Let us note that, V77 £ [0, 1[, we have 



sup \a j (u)-jn\=0(N(uj) a - 1 ), (7) 
je{i,...,N(u)<*} 



and 

lim sup \\<f>j(w,') - (t>j(-)\\H = 



with 



|sin(^x) if xe[0,d] 
if i>d. 



This result means that in the limit of large number of propagating modes the low order propagat- 
ing modes are very similar in shape to those of a perfect bounded waveguide with pressure-release 
boundary conditions at x = and x = d. This approximation does not hold anymore for high or- 
der propagating modes, but the results of Lemma [l . 1 1 mean that the distribution of solutions of (|6| 
is closed to the distribution of the eigenvalues of the transverse Laplacian associated to a perfect 
bounded waveguide with pressure-release boundary conditions at x = and x = d. 

Continuous part of the decomposition For 7 £ (— 00, k 2 (u>)), we have [21] 
7 (w, a;) = 

{A 1 {u))s\n{r]{ijj)x/d) if < x < d 

A 7 H(sin(r ? H)cos(eH^) + fMcos(^H)sin(eH^)) if d < x, 



where 



r}((jj) = d^n 2 k 2 (uj) - 7, f (w) = dy/k 2 (u}) - 7, 



and 



7t(£ 2 (cl>) sin 2 (?7(u;)) + r] 2 (oj) cos 2 (?7(w))) 



Let us note that </> 7 (w, .) does not belong to H so that (y, <^ 7 (w, .)) H is not defined in the classical 
way. In fact, we have 

(y,<£ 7 (w,.)) = lim / y(x)4>ry((j,x)dx 
where the limit holds on L 2 ( — 00, fc 2 (w)). Moreover, according to the following Plancherel equality 

N(ul) 2 f k2(hj) 2 

HvllH = l|n«(-oo,+oo)(i/)||^= \(v^M-))„\ 2 + / \(y,^,.)) H \ 2 dj, 

the map which assigns to every element of H the coefficients of its spectral decomposition 
9 W : H — > W 

is an isometry, from iJ onto = C N ( U > xi 2 (-oo, £; 2 (w)). 

2 Mode Coupling in Random Waveguides 

Before describing the time-reversal experiment in our randomly perturbed waveguide model, we need 
to understand how the wave is perturbed during the propagation through the medium. In this section, 
we study the random effects produced on the modal decomposition of p(tu,x,z) propagating in the 
perturbed section [Q,L/e]. 

Using the resolution of the identity 11^ associated to Pekeris operator R(uj), we have 

p(u,x,z) = Yp j {u3,z)^j{u,x)+ I p J (uj,z)(p 1 (uj,x)d'j, 
j=i J -°° 
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where p(u>, z) — Q u ,(p[u), ., z)) and 6^ is denned in Section 

For the sake of simplicity in the presentation of the forthcoming asymptotic analysis, we will 
restrict ourself to solutions of the form 

p{u),x,z) = ^2 Pj{u,z)<t>j{u>,x) + I p J (uj,z)cp 1 (uj,x)d'y. (8) 
j'=i ^ 

This assumption is tantamount to neglecting the role played by the evanescent modes during the prop- 
agation in the random medium. Nevertheless, as it has been observed in |15l 110] that these modes 
play no role in the refocusing process. The reason is that these modes only imply a mode-dependent 
and a frequency-dependent phase modulations without remove any energy from the propagating and 
radiating modes, However, the dispersion phenomena are compensated by the time-reversal mecha- 
nism. Moreover, we assume that e <C £ and therefore we have two distinct scales. We will consider in 
a first time the asymptotic e goes to and in a second time the asymptotic £ goes to 0. 



2.1 Coupled Mode Equations 

According to the pressure field decomposition pi), we give in this section the coupled mode equations, 
which describes the coupling mechanism between the amplitudes of the two kinds of modes, prop- 
agating and radiating modes. In the random section [0, L/e], p(u>,z) satisfies the following coupled 
equation in U% = C N ^ x £ 2 (£, k 2 {uo)). 

d 2 - N(UJ) 

— pjfa, z) + $MPi(«. z ) + V^H £ C%{z)pi(u, z) 
z i=i 



+ sfek 2 {u) I C^,(z)p Y (uj, z)di = 0, 



2 N( U ) 



(9) 



d 

p 7 (u),z) + jp^(u,z) + ^/ek 2 (uj) C^(z)pi(u; , z) 



dz 2 

i=i 



+ ^ek 2 (uj) C",(z)p Y (oj,z)dy = 0, 



where the coupling coefficients C w (z) are defined by: 



C ji( z ) = (</>3(u,.),<l>i(u,.)V(.,z)) II = / 0j(u,x)(j)i(uj,x)V(x 7 z)dx, 

Jo 

qi 1 {z) = C^{z) = (if> j (u,.),^(u,.)V{.,z)) H = [ <p j (oj,x)^(oj,x)V(x,z)dx, (10) 

Jo 

f d 

Cyy( z ) = (<P'y(u,-),<P 1 '( UJ ,-)V(.,z)) H = / <p 1 (uj,x)(j)y(uj,x)V(x,z)dx. 

Jo 

Next, we decompose the wave field p(u, z) using the amplitudes of the generalized right- and 
left-going modes a(w, z) and b(u, z), which are given by 

ftfaz) = —)=(a 3 (u,zy^ +b j (u,z)e- i ^ u >), 



d 

dz 



%[u, z) = ^7 1 / 4 (a 7 (w, z)e l ^ z ~ S 7 (cj, zje - *^*) 



Vj G {l, . . . , AT (a;)} and almost every 7 € (£, k 2 {w)). From ([9]), this decomposition allows us to obtain 
a first order differential system instead of a second order one, so that we obtain the coupled mode 
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equation in W£ x W£ for the amplitudes (a, b) , 



dz 
d 
dz 



a(ui,z) = \/fH a °(w,z) (a(u J ,z)) + ^eH. ab (u J ,z)(b(u J ,z)) 
b(uj, z) = v / eH ba (a;, z)(a(uj, z)) + </eH bb (oJ, z)(b(uj, z)) . 



(11) 
(12) 



This system is complemented with the boundary conditions 



a(oj, 0) = ap(o;) and b I u>, 



L 



where 



= ^^(^ ■),U-)) H e-^ )Ls , Vi e {1, • . ■ , N( U )}, 
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2 

1/4 



(13) 



-<$ e (w,-)> 7 (w)) ff e~^ Ls , for almost every 7 G (£,fc 2 (w)). 



For j G {l, ... , iV(w)}, 0^0(^0) represents the initial amplitude of the jth propagating mode, and 
for 7 G (£, fc 2 (w)), a 7i o(w) represents the initial amplitude of the 7th radiating mode at z — 0. The 
initial conditions for the right-going mode ao(wo) comes from ([3|, and the second condition for the 
left-going modes means that no wave is coming from the right homogeneous waveguide. The coupling 
operator H QQ (a;,z), H Qb (w,z), H &a (w,z), and H b6 (u;,z) in (ITT) and pi are defined by: 



Hj a (cj, z)(y) = Hf(uj, z)(y) = 



ik 2 (uj) 



N(u>) 



OK*) 



--Vie 



= <(^,z)( y ) = ^Mr V °^ z) — - 



N(u 



2/;e 



fc V) C",(z) 



7 l/4yl/4"^' 



2 v^HflH 



-i(ft(o;)+/3.,-(w))z 



[E 



fe2(w) <V(z) 



we 



-i(A(w)+ N /7)« 



s 



7 i/4yi/4' 



Let us remark that we have the following global conservation relations 

\\a(u,z)\\ 2 n „ - ||6(w,z)||^ = ||S(w,0)||^ - ||6(«,0)|& ? Vz e [0,L/e] 
||a(w,L/c) II 2 ,. + ||6(^,0)|| 2 = R W ,0)|| 2 



« i i 

However, in our context, the coupling mechanism between the right- and the left-going mode is not 
very convenient to study the time-reversal experiment. The asymptotic behavior of the whole coupling 
mechanism between the right- and the left-going in a random waveguide with a bounded cross section 
is carried out in [13] , but this study leads to technical difficulties because of the waveguide geometry 
in our context. Consequently, for the sake of simplicity we introduce in the following section the 
forward scattering approximation, which allows us to neglect the coupling mechanism between the 
right- and the left-going under certain conditions. 
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2.2 Propagator and Forward Scattering Approximation 



Before introducing this approximation, let us define the rescaled processes according to the size of the 
random section [0, L/e], 



a c (w, z) = a^uj, —J and 6 e (w, z) — b[ui, — ^ 



which satisfy in the rescaled coupled mode equation 

d^, , 1 



jj-qq 



-a e (w, z) 
-b e (uj,z) 



(V) (?Kz))4h« 1 ( W> ^) (&>,*)) 



1 



H 



6 a 



w,-) (a e (w,z)) + 



1 



H 



/,/, 



(& £ (c,z)), 



dz ' ' ' y/e 

with the two-point boundary conditions 

a e (uj,0) = a e (u) and b c {uj 1 L) = Q. 

The propagator P e (w, z) is defined has being the unique solution of the following differential 
equation 



dz 



P(u;,z) = (u, -) P e (w,z) with P e (w,0) = /d, 



so that 



a e (w, z) 
S e (w,z) 



P>,z) 



a e (w,oy 

b%io,0)_ 



According to the symmetry of H(w, z) the propagator has the following particular form 

P (u,z) - 



_P°(w,z) P°( w ,*) 

where, P°(w,z) and P'(w,z) are two operators representing respectively the coupling between the 
right-going modes and the coupling between the right-going and left-going modes. 

The forward scattering approximation is widely used in the literature. In this approximation the 
coupling between forward- and backward-propagating modes is assumed to be negligible compared to 
the coupling between the forward-propagating modes. The physical explanation of this approximation 
is as follows. The coupling between a right-going propagating mode and a left-going propagating mode 
involves coefficients of the form 



nC%(0)C%(z)] cos ((AM + Pi(u))z)dz, 



where the coefficients C u (z) are defined by (10 1, and the coupling between two right-going propagating 
modes or two left-going propagating modes involves coefficients of the form 



E[qi(o)C5J(*)]cos ((AM - PM) z )dz 



V(j, I) € {l, . . . , N(uj)} 2 . The forward scattering approximation consists in assuming that 



E[C$(0)C$(z)] cos (08,(w) + ^(u))x)dz = V(j,Z) e {l,...,N(u)} 2 , 

JO 

so that P*(w, z) — 0, that is there is no coupling between right-going and left-going propagating 
modes. Therefore, this approximation holds if the power spectral density of the process V, i.e. the 
Fourier transform of its z-autocorrelation function, possesses a cut-off wavenumber. We refer to 
|12L I15j for justifications on the validity of this approximation. As a result, under this approximation 
we can neglect the left-going propagating modes in the asymptotic e — >• 0, and then consider only the 
simplified coupled amplitude equation on [0, L] 



— a e (uj, z) 
dz 



(u>, (a e ([uj, z)) with a e (ui, 0) = a (uj). 
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Time-reversal 
mirror Vu 



Time-reversal 

mirror T>m 



(a) 



(b) 



Figure 3: Representation of the time-reversal experiment. In (a) we represent the first step of the 
experiment, and in (6) we represent the second step of the experiment. 



Finally, we introduce the transfer operator T^' e (o;,z), which is the solution of 
d 



dz y/e 



UJ, 



T«' e (w,z) with T 4 ' e (w,0) = Id. 



(14) 



From this equation, one can easily check that the transfer operator T^ £ (w, z) is unitary since H aa is 
skew-Hcrmitian and 

Vz > 0, a e (uj, z) = T 5 ' £ (w, z)(c%(u})). 



3 Time Reversal in a Waveguide 

Time-reversal experiments with sonar in shallow water [181 I23| were carried out by William Kuperman 
and his group in San Diego. This experiment is carried out in two steps. In the first step (see Figure 
[3] (a)), a source sends a pulse into the medium. The wave propagates and is recorded by a device 
called a time-reversal mirror. A time-reversal mirror is a device that can receive a signal, record it, 
and resend it time-reversed into the medium. In other words, what is recorded first is send out last. 
In the second step (see Figure [3] (6)), the wave emitted by the time-reversal mirror has the property 
of refocusing near the original source location, and it has been observed that random inhomogencitics 
enhance refocusing [3J [5J 03 [HI [HI EH [H] . This experiment has already been analyzed in waveguides 
with bounded cross-section in (TUl Chapter 20] and [HI HI]. However, in contrast with all these 
results where the random medium improve the refocusing, we show in this section in a context of 
underwater acoustics that the random inhomogencitics deteriorate the refocusing property. In [T] the 
authors observe such deterioration but in their context it is induced by viscosity effect in the wave 
propagation model. We show in this section that this effect is simply induced by the inhomogeneities of 
the propagation medium through the coupling mechanism between the propagating and the radiating 
modes. 



3.1 First Step of the Experiment 

In the first step of the experiment (see Figure |4| , a source sends a pulse into the medium, the wave 
propagates and is recorded by the time-reversal mirror located in the plane z — L/e. We assume 
that the time-reversal mirror occupies the transverse subdomain T>m C [0, d] and in the first step of 
the experiment the time-reversal mirror plays the role of a receiving array. The transmitted wave is 
recorded for a time interval 7*] and is re-emitted time-reversed into the waveguide toward the 
source. We have chosen such a time window because it is of the order of the total travel time of the 
section [0, L/e]. 

In this paper, the source profile ~$! e (t, x) is given, in the frequency domain, by 



« v ' ' e 9 V e« 



V (j)j(u),X Q )<l)j{u,x) + / 

J( 



(-s-£)u(«,fc 2 (^)) 



4> 1 (uj,x )(j) 1 (ui,x)d'y 



(15) 
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Figure 4: Representation of the first step of the time-reversal experiment. 



with q > 0. The restriction q > allows us to freeze the number of propagating and radiating modes, 
introduced below, and gives simpler expressions of the transmitted field. Let us note that S can be 



arbitrarily large and £ can be arbitrarily small, so that the transverse profile ( 15 1 is an approximation 



of a Dirac distribution at Xq, which models a point source at xq. Moreover, ^/( "~"° ) is the Fourier 
transform of f{t q t)e~ lu}ot , which is a pulse with bandwidth of order e q and carrier frequency U)q. In 
this paper, we are interested by a source emitting a broadband pulse, that is for q £ (0,1). For 
our broadband source term with pulse width of order 1/e 9 , smaller than the propagation distance, 
the propagating modes are separated in time by the modal dispersion. With this kind of source, 
waiting long enough to record all the train of pulses at the time-reversal mirror, one can observe a 
self-averaging effect on the refocused pulse. This statistical stability implies that the refocused pulse 
does not depend on the particular realization of the random medium. This phenomenon has been 
widely studied in different contexts and there are many references about it [3J HJ |5J [TOl [H| ■ The 
case q = 1, that we do not treat in this paper corresponds to the narrowband case. In this case the 
order of the pulse width is comparable to the propagation distance. Consequently, the modes overlap 
during the propagation and then the statistical stability of the time-reversal experiment depends on 
the number of propagating modes jTUl [H] • However, for the sake of simplicity, we consider only the 
case 9=1/2 but the following analysis can be carried out Vg € (0, 1). 

According to [TO] HH [H>] , the evanescent part of the wave field decreases exponentially fast with the 
propagation distance. For more convenient manipulations in the study of the time-reversal experiment 
we assume that the source location L$ is sufficiently far away from so that the evanescent modes 



generated by the source are negligible. With this assumption and using (13), we can assume that the 
incident pulse coming from the left is given, at z = 0, by: 



*4c(*,z,0) = 



1 

27 



<t>j{u,x) 



k 2 (u) 



l 7,0 



,1/4 



t du, 



where 



Vje {l,...,N(u>)}, 

«7,o( W ) 



1/4 



2ei 



(16) 



(17) 



for almost every 7 6 (£, k 2 (uj)). Let us remark that this assumption is not restrictive and all the results 
of this paper are valid for any L5 < 0. Indeed, according to Proposition 4.2 in [TS], in the asymptotic 
e — > 0, the information about the evanescent part of the source profile are lost during the propagation 
in the random section [0, L/e], and therefore they play no role in the pulse propagation and in the 
time-reversal experiment. An efficient way to do not loss the information about the evanescent part 
of the source term has been developed in [5D] and studied in [2] . 
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Figure 5: Representation of the second step of the time-reversal experiment. 



Finally, according to Section |2j the wave recorded by the time-reversal mirror is given by 



Ptr \t,x, 



f 



OJ — LJq 



k 2 . 



l 



^T^' e ( W ,L)(a( W ))0 7 ( W ,x) e V7t d7e 



/« 7 

where T^' e (u>,L) is the transfer operator solution of (14) 



3.2 Second Step of the Experiment 

In the second step of the experiment (see Figure |5| , the time-reversal mirror plays the role of a source 
array, and the time-reversed signal is transmitted back. Now, the source term is given by 



with 



F e TR (t, x, z) = -fx R {t, x)8(z - L/e)e z , 
f TR (t,x)=pl (| - t,x, ^ G^h - et)G 2 (x), 



where 

= l[ta,td(*) and G 2 {x) = 1 Vm (x). 
Here, G\ represents the time window in which the transmitted wave is recorded, and G 2 represents 
the spatial window in which the transmitted wave is recorded. In our study, we are interested in the 
spatial effects of the refocusing, so we assume that we record the field for all time at the time-reversal 
mirror, that is the source has the form 



fj- R {t, x) = p tr - t, X, G 2 (x). 



Now, we are interested in the propagation from z = L/e to z — 0. The decomposition with respect to 
the resolution of the identity n w associated to R(u) (see Section |1.1| gives 

p TR (u,x,z) = ^pfile-^^^x) + / M^ e -^ >7x)d7j 



=1 V/^mM 



7 



with 



b m (u,L) 
b 7 (uj,L) 



vV4 



</Ti?(w,-)^m(w,.)> ff , 

(/rfi( w v),^( w ^)) H 
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in W£ . Then, at the source location z = L s , we obtain 



with 



&(w,0) = (T^)*(uj,L)(b(uj,L)), 
where (T*' e )*(cj, z) stands for the adjoint operator of T^' e (w, z). Consequently, one can write 

PTfl(^x,L s ) = (T^(w,L)(6 x (w)),6(w,L)) HU ,, 



where 



and 



6*,nM= -^==M^x)e- i ^ La and 6 Xi7 (w) = ^^(w.^e"*^, (18) 



4\A 

in H| , where A e (w) is defined by 



6>,L) = J-/f^)e^(^( W ,I)(a( W )),A;( W )) HB 



A» 



and with 



1/4 

M rs (cj) = / G2(x)</v(cj, x)(j) s (tj, x)dx 
Jo 



(19) 



for (r, s) e ({1, . . . , N(u>)} U (£,fc 2 (cj))) . (M rs (w)) represents the coupling produced by the time- 
reversal mirror between the modes during the two steps of the time-reversal experiment. Consequently, 



ptr^x iLs ) = ^/(^)^ 1 (u^(^i)(«(^).^H) I A £ H) J1 „ Mr 

Here, we consider the tensorial space W£ <g> — |a <g> fi, (A, fi) G (%^ ) 2 j, with (A ® ^) r;s = A r ^ s 

for (r, s) e ({1, . . . , JV(w)} U (£, fc 2 M)) and V(A, /x) e (?^) 2 . This space is equipped with the inner 
product defined by 

JV(«) JV(w) „ fe 2 



( x >v)n«®n? = 12 x oiWi+ 12 / ^-i'WsH 
5 4 i,J=i j=i J? 

/■/s 2 (w) "(") /-*: 2 (w) /-*: 2 (w) 



/« ,=i ^ J ? 

V(A,/x) e ("Hj Finally, the time-reversal kernel U e ' e (w, L) is defined by 



U«> £ ( W ,L)( 2/ 1 , 2 / 2 ) = T^( W ,L)(yi)®T^( W ,i)(y 2 ) (20) 
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V(y 1 ,y 2 ) G {Ti^) 2 , describing the two steps of the time-reversal experiment through the random 
medium thanks to the transfer operator T^' e (w,L) satisfying (14). 



We study the refocused wave in a time window of order l/y/e comparable to the pulse width, 
and centered at time t t, s /e, which is of the order the total travel time for a distance of order 1/e. 
Consequently, we will study the refocusing of the refocused wave at the source location z = L$ given 
by 



Vtr 



^obs 



1 



8^ 



(U«' e (w, L)(o(w), 6,M), A e M)^^e 



(21) 



duj. 



where a{oj) is defined by (16 1 and (17), b x (w) is defined by (18), and A £ (w) is defined by (19 1 
In what follows, we consider a time-reversal mirror of the form T>m = [^1,^2] with 



dn = d 



m + Kc'd2 and d 1 = d M - d u 



where dp,{ G (0,g?), (0^2,^1) G (0,+oo) 2 , and G [0,1]. Here, A oc = 27rc/(niWo) is the carrier 
wavelength in the ocean section [0, d] of the waveguide. The time-reversal coupling matrix are therefore 
given by 



Mfj(w) = (d 2 - d^Aji^Aiiu) 



cos [ (cr^w) - j Sine ^(0-3(0;) - <T l( UJ )) d ' 2 2d dl 

cos [ (o-j(w) + <J i( UJ )) d ' 2 2^ I sinc ( fo'M + ^/M)^ 



(22) 



2d 



for (j, Z) G {1, . . . , N(uj)} 2 , where Aj(uj) and Oj(w) are defined in Section [l.l| We give only the 
coefficients MjiiuS) for G {1, . . . , N(u>)} 2 , because in what follows only these terms will play a 
role. The parameter olm represents the order of the magnitude of the size of the time-reversal mirror 
with respect to the wavelength in the ocean cross-section [0,d]. In fact, we will see that the size of 
the mirror plays a role in the homogeneous case only when it is of the order the carrier wavelength 
A oc = 2nc/(niUJo)- 

In the following section we study the transverse profile of the refocused wave in the continuum 
limit iV(wo) 3> 1 of a large number of propagating modes, which corresponds to the regime luq +00. 
However, we know that the main focal spot must be of order A oc , which tends to in this continuum 
limit iV(wo) S> 1. Consequently, we study the transverse profile of the refocused wave in a spatial 
window of size A oc centered around xq. 



3.3 Refocused Field in a Homogeneous Waveguide 

To understand what are the effects produced by the random perturbations of the propagation medium 
on the time-reversal experiment, we study first the refocused wave obtained in a homogeneous waveg- 
uide. Let us consider the refocused wave in a time window of order \/-\fe, which is comparable to 
the pulse width, and centered at time t b s /e, which is of the order the total travel time for a distance 
of order 1/e. In this section we assume that the medium is homogeneous, so that T^ e (u,L) = Id. 
Then, the refocused wave at the original source location is given by 



Vtr 



V e Je J 4 ^— ' 

j,m—l 

x 0j( w o, x o )0 m (w o , x)K J ^ L * / ^ J — j= t J 



where 



^C.lM = K Tl(")KZl{") = e^'M-^o^lr, (23) 
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and L — S . Consequently, in the asymptotic e — > 0, we can observe a refocused wave only for a 
finite set of times given by 

t jm = t 1 + (p m (u )-P' J (u ))L. (24) 

For m 7^ j, we obtain 

V e ye / 

+ 0(Ve). 

At time tj m (j ^ m) one can observe only the mth mode, emitted by the time-reversal mirror during 
the second step of the experiment, coupled with the jth modes recorded by the time-reversal mirror 
during the first step. This coupling is produced by the time-reversal mechanism through the time- 
reversal mirror and characterized by the coupling matrix Mj m {u)o). Moreover, let us note that the 
refocused wave shape is dispersed by the kernel K"° L (t) during the first step of the experiment and 
by K^ L (—t) during the second step. 
Now, for t i, s = t\ we obtain 

pm^j + ^,x,L s ) = e~ iu °^ f{-t)H^{^x) + 0{yfe), 

where 

H x"{uo,x) = - M w -(w o )0j(w o ,a;o)0 3 (wo,a;). 

Here, we have a contribution of all the propagating modes. The refocused wave is a superposition of 
modes where each mode is coupled with itself by the time-reversal mirror through the terms Mjj(u>o). 
We describe in the following proposition the transverse profile of the time-reversed pulse in a very 
simple way in the regime ujq +oo, which correspond to the continuum limit N(ujq) 3> 1 of a large 
number of propagating modes. 

Proposition 3.1 For olm £ [0?1); the transverse profile of the refocused wave in the continuum limit 
N(ujq) 3> 1 is given by 

The width of the focal spot is therefore given by the diffraction limit X oc /(29). 

The proof of Proposition |3 . 1 1 is given in Section |4~T| As a result the sine function describes the asymp- 
totic transverse profile of the refocused wave in the continuum limit. 

The sine profile has already been obtained in different contexts in time reversal to describe the 
refocused transverse profile. Moreover, the size of the focal spot is of order A oc /(26>), where A oc is 
the carrier wavelength of the ocean section of the waveguide, and 9 depends on the contrast of the 
refractive index between the ocean section and the bottom of the waveguide. In [5) [TU1 [HI [H] for 
instance the random perturbations of the medium improve the time-reversal refocusing, but we will 
see in what follows that this statement is no more true in our context if the mode coupling mechanism 
between the propagating and the radiating mode is not negligible. 



3.4 Limit Theorem 



To describe the effects of the random medium on the time-reversed wave (21 1 we need to know the 
asymptotic distribution of the process U^' € (w, .), defined by ( pp| as e goes to and £ goes to 0. First, 
let us remark that V(y\ y 2 ) € (H^) 2 , with H<£ = C N{uj) x £ 2 (£, k 2 {uo)), 



\\V^(^z)(y\y 2 )\\ 2 



Wz > 0, 
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and then let us introduced some notations. Let r„ 



ly 1 ® y 2 \\n^®H^ 



B ry ,u«®u t = {\&U%®iq,\\\\\u t ®u t <r y ) 



the closed ball with radius r y , and {g n ,n > 1} a dense subset of B ryt -H^®W£- We equip B ry ^®W£ 
with the distance d$ r u u, 9H ^ defined by 



+00 ^ 

^E^ ( A --",5n)^ 



V(A,/x) e {Br y ,H^®H^) , so that (B-upd^^^) is a compact metric space. 

In the following theorem, we give only the drifts of the infinitesimal generators because only this 
part is of interest in what follows. 

Theorem 3.2 V(y 1 ,y 2 ) € {%^) 2 , the stochastic process U^' £ (w, ^(y 1 ,!/ 2 ) converges in distribution on 

C([0,+oo),(B ry ^^^,di3 r w <v 8W u<)) as e — > to a limit denoted by V^^uj,.)^ 1 , y 2 ), unique solution 

0/ a well-posed martingale problem on ® "H^ starting from y 1 y 2 . Moreover, V(y 1 , y 2 ) G (Hq) 2 , 

the stochastic process U^(w, -Xy 1 , y 2 ) converge in distribution on C([0, +00), (B ry ^(g,-H^ , dg^ )) 

as £ — » to a fe'rm'i denoted by U°(w, .)(y\ y 2 ). T7ws fe'rm't is £/ie unique solution of the well-posed 
martingale problem on Hq ® "Hq starting from y 1 ® y 2 , and w'tft drz/i yzuen 



where 



N(u>) 



iV( w ) 

+ 2 E R» + r »H - ( r }» + r «H - 2f »)] (Pudu,, + U 3l d^ 

+ 2 E / t r , c » - T )M ( U ^ d u„ 2 + u^djj-y^ 

3 = 1 Ji 

1 ,fc 2 ( w ) W M _ 
+ 2 / E t r »M - r »H] + ^%^7i 

v C / 1 



N(u>) 



,fc 2 ( w ) 

E / r ij( w )( C/ i72^ 2 -Uj-ndfj— )d l2 

fe 2 (^) W M __ 

I] rf ; ( w )(c/ 7i; a^ ii - c^i%^)d7i, 

z=i 



1G 



and 



3,2=1 

\ E [Af H - A»] (l^fy, - t^) 

3,2=1 

1 ,fc 2 ( w ) 

- £ / [A» " *A 3 »] + [AJ( W ) + iA»] 



i=i ^ 
fc 2 (w) n(u) 

E [ A *H + iAfM]^!^, + [AfM - iAf(w)]t^r%^. 

Z=l 



The proof of Theorem 3.2 uses exactly the same techniques as the ones developed in [T^l [TB]. They 
are based on the perturbed-test-function method introduced in [19 and martingale techniques. Here, 
we have considered the complex derivative with the following notations. If t/ = U 1 +iU 2 € Hq ®%o, 
we have (Z7 1 , U 2 ) G ® £o ) 2 , where ^ = R w (") x i 2 (0, A: 2 (o;)). Then, the operators % = (du r ,J 
and djj = (djj — ) are defined by 



1 



1 



d u = 2^ 9ul ~ idu ^ and d u = 2^ 9ul + i9 ' 
with V/ G C 1 ((^o ® £o ) 2 > R ) and VA = (A 1 , A 2 ) G ® ££) 2 

E [ E + E / ^ 2 d u? j(v\v 2 )d l2 

n=l,2 3,i=l 3=1 J £ 



'r/ 2 J 



fe 2 (u) pk 2 (u>) 



E <A",^ / ( U 1 ,, 2 )) e ^ e ^^/(, 1 , v 2 )(A) 



which is the differential of /. Moreover, r c (ijj), r s (w), T (w), A c (ijj), and A s (w) are defined as follows: 
V(7,Z)e{l,...,JV(w)} 2 andj^i 



fc 4 (w) 



E[qj(o)qi(*)] cos ((AM - ftM)*)**, 



;=i 

'#3 



r?z(w) 



fc 4 (w) 



+ 00 



2/3 i (u>)/3 ; (w) J 

N{lu) 

r 33H = -E r 32H' 



E[qj(0)C*i(js)] sin - Pj(u>))z)dz, 



i=i 



and V(i,0 G {l,...,iVH} 2 



A» 



2&(w)#(w) Jo 



fc 2 (^) fc 4( w ) r+oo 



E[C^,(0)C^,(z)] cos ((VY - p j {u)))z)dzeh , > 
E[C^,(0)C^,(z)] sin ((VV - ftM)*)***/, 
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where the coefficients C w (z) are defined by (10 1. 

From Theorems |3.2| we have the following proposition about the autocorrelation function of the 
transfer operator for the two steps of the time-reversal experiment. 



Proposition 3.3 V(y ,y ) € Hq xHq anrfVA S x^jj, i/ie autocorrelation function of the transfer 
operator for the two steps of the time-reversal experiment as e — >• and £ — >■ is given &y 



lim lim E 



<U«<>, Ljfc, 1 , y 2 ), \) HU9H „] = E [<U>, L)(y\ y 2 ), A) 



j,l=l j,m=l 



1 (")-rL m (")-A^("))i+2( r ™ m (")- A ™M) i ^T ?/ 2 a — d-y 



3=1 



70 m=l 



fe 2 (aj) pk 2 (u>) 



yiyLX^yd-ydj'. 



Here, 

Qjmi") =\ [q» + - (T*» + T l mm {u) - 2r},(w)) - (A» + A^(w))] 

+ 1 [r^ m («) - l*» - (A^( W ) - a?(«))] . 

w'i/i Tj(u>, z) is the solution of the coupled power equations 



dz 



(25) 



and Tj(oj, 0) = <5jj. 



Let us note that, 



Tj(u Q ,L) 



lim lim E 



Tf'^WB.L)^)! 1 



(26) 



with y l j — Sji and y*, = for 7 6 (0, fc 2 (w)), is the asymptotic mean mode power of the jth propagating 
mode of the transfer operators at distance z — L. The initial condition y l means that an impulse 
equal to one charges only the ?th propagating mode at z = 0. Equation (25 1 describes the transfer 



of energy between the propagating and the radiating modes through the energy transport matrix 
T c (uj) the dissipation coefficients A c (ui). These dissipation coefficients resulting from the coupling 
between the propagating and the radiating modes are responsible to the radiative loss of energy of the 
propagating modes into the ocean bottom. As we will see in the following section the deterioration 
of the refocusing is due to these radiative losses which are caused by the random perturbations of the 
propagation medium. 



3.5 Refocusing in a Random Waveguide 

Using the asymptotic analysis developed in the previous section, now we are able to describe the 
asymptotic mean refocused wave. However, let us note that we only need to know the asymptotic 
mean refocused wave since the refocused wave is self-averaging, which means the refocused wave 
converges in probability to its asymptotic mean value. We refer to |161 for instance for a complete 
proof of the self-averaging property. Let us note that the self-averaging property of the time-reversal 
process has already been observed in many context [3J [5] EH ED El]- 
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Using the change of variable u) = loq + \feh, the refocused wave is given by 



Ptr 



(t 



--,x,L s )e 



1 

8^ 



f(h)e' 



tl-tpbs 



(U^ e (u; + Veft, L)(o(wD + </eh),b x (wo + V~eh)), A e (w + V^h)) nU0+Vi H 



Ah. 



Using Proposition |3.3| we obtain 



E 



—j=,x,Ls 



1 [rt^ 7 

£ / Pm\UQ) i(p n 

j,m—l V J 



( W o)-ft(^o)) J 



x M m j(uj )K^° n L * f 



i E 



U f m ( w o , i ) (a ( w ) , b x ( w ) ) 



(27) 



where ^ are defined by ( 23 1 . Let us note that there is no radiating part in the expression of the 



refocused wave because all the radiating components of the refocused wave involve a term of the form 



uniformly bounded in x and y on bounded subset of [0,+co) 2 . Moreover, we cannot observe the 
recompression of the radiating components by time reversal because it holds only on a set with null 
Lebesgue measure. 



From ( 27 1 , we obtain for m ^ j 
>t 



lim lim E 



Ptr 



e 



t 



, x, L f 



-i(/3 m (u l0 )-f3 j (a, )){-L s + i 



OJQ 

j,m. 



L * /(-*) ' ^j(^o,xo)(f> m (^o,x), 



where the times tj m are defined by (24 1. Then, at all these times we can observe the shape of the 



refocused waves obtained in an homogeneous medium, but with the damping terms e ( ^" l ( w °) i . As a 
result, the amplitude of the coherent refocused waves at times tj m decays exponentially with respect 
to the propagation distance L, and therefore becomes negligible for long propagation distance. More 
precisely, as we will see in what follows, even for t f, s = t\ the amplitude of the refocused wave will 
decay exponentially fast, but in this case the decay rate is smaller that in the case t b s — tj m (m =/= j). 
Now, for t b a = ti, we have from Proposition |3.3| a contribution of all the propagating modes 



lim lim E 

J-i-Oe^O 



Ptr 



-.,x,L< 



= /(-*)■ 4 Yl M v{ U) o)T-{ujo, L)(f>i(uj ,x Q )(f>i(uJo, x), (28) 

3,1=1 



where T\ (ujq , L) are the asymptotic mean mode powers ( 26 ) satisfying the coupled power equations 



(25 1, and Mjj{uo) is defined by (22 1. From this last expression one can see that the refocusing of 



time-reversed wave is closely related to the mean mode powers propagation described through the 



coupled power equations ( 25 ) . We will sec in the following sections that the effects of the radiative 



losses into the ocean bottom, and caused by the random inhomogcncitics, affect the refocused wave 
in two ways. First, the amplitude of the refocused wave decay exponentially with the propagation 
distance. This exponential decay of the propagating mode powers is described in |17M15j . The second 
effect is the loss of resolution of the transverse profile. Let us recall that for a waveguide model with 
a bounded cross section the energy is preserved, so that one can observe the classical time-reversal 
superresolution effect mainly through a side-lobe suppression produced by the time-reversal mirror 

nana. 



3.6 Exponential Decay of the Refocused Wave Amplitude 



The exponential decay rate of the asymptotic propagating mean mode power ( |26| is carried out in 
|15j . but let us recall the result after introducing some notations. 



19 



Let us consider 



S 



JV(i 



Jo) = R Nt -^'\ X.j>Q Vj G {1, . . . , iV(w )} and ||X|| 



21™("o) 



with (x, y) 



Efli X 3 Y 3 for e (^ ( " o) ) 2 , and 



,A 



)(^o))- 



(29) 



Theorem 3.4 Le£ us assume that the energy transport matrix r c (wo) is irreducible. Then, we have 

1 



lim — In 



4=1 



-Aoo(a;o) 



with 



A oo ( W0 )= inf ((-r c ( Wo )+A^M)X,X) 



XG5 



+ 



>0. 



Unfortunately, A oc (o;o) is not easy to compute, but however we have the following inequalities 

1 



< A r , 



je{i,...,N(u )} 



A-(w ) < A 00 (wo) < A(w ) 



N(u ) 



(30) 



Let us give two examples for which A oc (o;o) can be computed explicitly. First, if we assume that 
the energy transport matrix T c (uj ) can be replaced by ^T c (uj ) with r«l, that is we assume that 
the mode coupling is strong. In this case one can show [15j . using a probabilistic representation of 
Tj (ujQ, L) in terms of a jump Markov process, that 



lin i A L( w o) = A(w ), 



and 



lim L) = cxp ( - A(cj ) 



L 



The idea is that we have a strong mixing so that the decays rate averages. On the other hand, if we 
assume that the energy transport matrix r c (wo) can be replaced by tT c (u>q) with r< 1, that is we 
have a weak mode coupling. In this case, one can show |15j that 

lim A^(w ) = A min (uj ), 

The idea is that the modes coupling is too weak to provide any mixing so that the decay rate for the 
jth mode is given by A^(ujq). 

As we will see in what follows, the loss of resolution in the time-reversal experiment does not 
depend only on the mode coupling between the radiative and the propagating modes, but also on the 
energy transfer between the propagating modes, and described by the transfer matrix T c (u>o). 

3.7 Refocused Transverse Profile without Radiative Losses 

Before, studying the effect of the radiative losses on the focusing quality, let us recall the basic result 
in the case of negligible coupling between the propagating and the radiating modes |101 112j . 
In this section and the following one let us note 



N{u ) 



H%™{uq,x,L) = - ^2 M jj (u )T : l(uJo,L)(j>i(uj Q ,x )(pi(uJo,x), 



(31) 



3,1=1 



the transverse profile of (28 1, where Mjj(u>o) is defined by (22 1, and cum represents the order of 
magnitude of the time reversal mirror. 

In this section, we assume that the radiative losses in the ocean bottom caused by the random 
perturbations of the propagation medium are negligible, that is A c (wo) = 0. Consequently, the asymp- 
totic mean mode powers satisfy (25 1 without the radiating coefficient A c (wq). Moreover, according 
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to |1(J1 112] , using a probabilistic representation of the propagating mean mode powers in terms of a 
jump Markov, we have 

lim 7~'(w ,L) = T77— r- 

L-5.+00 J l\(U}o) 

This result describes the asymptotic behavior of the propagating mode powers for long propagation 
distances. In this asymptotic the energy is uniformly distributed over all the propagating modes and 
there is no loss of energy, in other words there is no loss of information. As a result, the refocused 



transverse profile (31 1 is given by 



JV(w ) ^ Af(w ) 

l^^o'^o^^L) = - M iAu )— ^ <f>i{u Q ,x )(j)i{u} Q ,x), (32) 
which is closed, up to a multiplicative constant, to the sine profile according to Proposition |3.1| 



However, let us note that in contrast to Section 3.3 for a homogeneous waveguide the resolution does 



not depend on the time-reversal mirror. Let us also note that this transverse profile corresponds to the 
projection of 5{x — xq) over the propagating modes, which corresponds to a point source localized at 
x = xq. In absence of energy loss, or information loss, the refocused transverse profile, which is closed 
to the sine profile, is then the best profile which can be expected for the time-reversal experiment in 
a random waveguide. However, to observe better refocusing properties, the waveguide setup has to 
be modified as described in [H] . 

3.8 Example of a Refocused Transverse Profile without Loss of resolution 

In this section, we show that the loss of resolution is not only due to the mode coupling between 
the radiative and the propagating modes. In fact, in the first example introduced in Section |3.6|thc 



resolution quality is not affected, we still obtain the profile (32) but with a damping term decaying 
exponentially fast, and with a decay rate depending on the mode coupling between the radiative and 
the propagating modes. 

If the propagating mode coupling process is stronger than the radiative losses, that is we consider 
^T c (ojo) with r C 1, we have 



r^o 7 ^^ ' L ) = tv^J) exp ( ~ A ( w o)^) , 



where A(wo) is defined in (30 1. Consequently, the refocused transverse profile (31 1 is given by 



N(u> ) ^ AT(w ) 

i ^™ 00 H x M (^o,^L) = -exp (- A(u) )L) Y n(ui ) ^ x )<j)i(u) ,x), 

-> 00 ._ x (u ) (=1 



which is the refocused transverse profile obtained without any radiative losses ( 32 1 but with a damping 



term given by the average radiative losses. Because of the strong mixing property the energy carried by 



the propagating mode is uniformly distributed as discussed in Section (3.6 1, but decaying exponentially 



fast, so that the transverse profile is the projection of S(x — xq) over the propagating modes. 

3.9 Loss of Resolution for the Refocused Transverse Profile in the Con- 
tinuum Limit 

This section is devoted to the study of the loss of resolution of the refocused transverse profile. 



Considering, a nearest neighbor coupling mechanism in ( 25 1 , we use the continuum approximation 



(N(u>o) S> 1) of this equation to describe the refocused transverse profile. 
3.9.1 Nearest Neighbor Mode Coupling 

In the following sections we investigate the loss of resolution of the refocused wave. To lighten the 



loss of resolution and in order to give simple representation of refocused the transverse profile (31 1, 



21 








z 







^ 






/ Ocean 




d 






Bottom 



Figure 6: Illustration of negligible radiation losses in the shallow- water random waveguide model. 

we consider the simplified coupled power equation 
d 



,T n {uq, z) = -A c n (uj )Tn(u} , z) + T c n ^ 1n (uj ) (7w_ 1 (w 0j z) - T^(ujq, z)) , 

~Tj(u , z) = r5_ Xi (wo) (Tj.^o, z) - 7> , z)) 

+ T c j+lj (uj ) {Tj +1 (u Q ,z) - TJ(oj q , z)) for ./ {2 V L}, 

£tK"o,z) = r§i(wo) (75(wo,«)-7?(wo,z)) , 



(33) 



with TJ(ujo,0) = Sji, that is we only consider a nearest neighbor mode coupling. As a result, only 
the highest propagating mode can still be coupled with the radiating modes since it is the closest 
mode to the continuous spectrum of the Pekeris operator. This simplified version of the coupled 
power equation can be rigorously derived using a band-limiting idealization |15l 117] , that is the power 
spectral density of the random perturbations in the transverse direction has a compact support. 
This simplification will be the starting point to study the refocusing properties using the continuous 
diffusion approximations of ( 33 ) introduced in [15] . These continuum approximations allow to exhibit 
in a very simple way the effects of the radiative losses on the transverse wave refocusing. First, let 
us recall what happened if the radiation losses are negligible, that is we neglect the coupling between 
the propagating and the radiating modes. Afterward, we describe the influence of the radiative losses 
in the ocean bottom on the time-reversal refocusing. 



3.9.2 Refocused Wave in the Continuum Limit with Negligible Radiation Losses 

In this section, we study the transverse profile of the refocused wave in the case where the radiation 
losses are negligible, that is A^r(u;) = 0. In this case, the two following propositions describe the 



refocusing properties of (31 1 



According to Proposition |3.1| let us recall that the size of the focal spot is of order the carrier 
wavelength of the ocean A oc = 2ttc/ (b^q), which tends to in this continuum limit N(loq) 3> 1. Let 
us also recall that the continuum limit N(uj q ) ^S> 1 corresponds to the case ujq /~ +oo. Consequently, 
we study the transverse profile (31 ) of the refocused wave in a window of size A oc centered around xq. 



Proposition 3.5 For olm £ [0jl); with negligible radiation losses, the transverse profile of the refo- 
cused wave in the continuum limit N(ujq) 3> 1 is given by 



lim 



_ IJOtM 

X 



WO, £o 



L\ = 



di 



sinc(27rx). 



The transverse profile of the refocused wave is studied using the continuous diffusion approximation 



model of the coupled power equation (33 1 introduced in [151 Theorem 6.5]. Without radiative losses 
the energy is conserved (see Figure [6| so that we obtain the best transverse profile, which is the sine 
profile. In the same way, we have the following result for a m = 1, that is the carrier wavelength 
and the magnitude of the time-reversal mirror are of same order. In order to study this case, let us 
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introduce some notations. Let £ = Um>i ^Mi where 

(oy)j G R M \ , and ^(x) = y^sin (j^) Vx G [0, d],Vj > 1. 



-M 




Let us remark that (<f)j)j is a basis of L 2 (0, d). 

Proposition 3.6 For olm — 1, £/ie continuum limit N(ujq) 3> 1, we have 



lim J? x 1 o (^,.,i)-^ o (w ,.,i)=0 



in £' , and where 

lim if* fcjQ) xq H — ^pic, = 9 2 — -sinc(27ra). 



These results are consistent with the ones obtained Section 3.7 and in [TUJ Chapter 20] and [12], 
where the authors have obtained the sine function for transverse profile. The most important fact is 
that, according to Proposition |3 . 1 1 the transverse profile of the refocused wave does not depend on the 
magnitude of the time-reversal mirror, we obtain the sine profile Vqm € [0, 1]. 

However, when the radiation losses are not negligible anymore the refocusing properties will be 
affected in two ways. First, the amplitude of the wave decrease exponentially fast with the size of the 
random section, and second, the refocusing quality itself is deteriorated by its losses. 

3.9.3 Refocused Wave in the Continuum Limit with Radiation Losses 

In addition to the exponential decay of the refocused wave amplitude, the second main effect of the 
radiative losses is a deterioration of the focusing quality described in the two following theorems. 



To study the refocused transversed profile (31 1 we use the continuous diffusion approximation model 



introduced in [T^l Theorem 6.3] of the coupled power equation (33). According to Proposition 3.1 let 
us recall one more time that the size of the focal spot is of order the carrier wavelength of the ocean 
A oc , which tends to in this continuum limit N(ujq) 3> 1 (<Jq +oo). Consequently, we study the 



transverse profile (31 1 of the refocused wave in a window of size A oc centered around Xq. 

Theorem 3.7 For oim € [0,1), the transverse profile of the refocused wave in the continuum limit 
N(u>q) ^> 1 is given by 



wht 



H(x,L)= / 7i(i, u) cos(2irux)du, 
Jo 

and 7l(L,u) is the solution of 



(34) 



l z TiM= L( a ^iL r ^ { u] 



with the boundary conditions: 

8 



Vz > 0. Here, 



. 71(z,Q)=Q, 7i(-M) = and 71(0,10=1, 
ou 



1 - 1 - ^ (On) 2 



a 2 d 2 

with a = 2 dfd°e 2 > ® ~ \A — V n !> = Jo Jo 7o(^i, ^2) cos Xi) cos (^x 2 )dxidx 2 . n\ is the index 
of refraction in the ocean section [0,d\, 1/a = l ZiX is the correlation length of the random inhomo- 
geneities in the longitudinal direction, and 70 is the covariance function of the random inhomogeneities 
in the transverse direction. 
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Figure 7: Illustration of the radiative loss in the shallow-water random waveguide model. 



The proof of Theorem |3.7| is given in Section [4. 3| Consequently, the transverse profile of the refocused 
wave can be expressed in terms of the diffusive continuous model, with a reflecting boundary condition 
at u = (the top of the waveguide) and an absorbing boundary condition at u = 1 (the bottom of 
the waveguide) which represents the radiative loss (see Figure [7]). 

Thanks to Theorem |3.7| we can now give a simple explanation of this refocusing degradation. 
First, let us recall that the quality of the refocusing depends much more on the high propagating 
modes that the lower ones. As illustrated in Figure [9] the absorbing boundary condition at u = 1, 
describing the radiation losses, degrades the information carried by the high propagating modes and 
it is getting worst as the size of the random section becomes large. As a result, the function 7i(L, u) 



in (34 1 plays the role of a low pass filter, which therefore degrades the refocusing quality as illustrated 
in Section 3.9.4 (Figure [9] and Figure 10). Consequently, the radiation losses degrade the quality 
of the refocusing in two ways: the amplitude of the refocused wave decays exponentially with the 
propagation distance, and the width of the focal spot increases and converges to an asymptotic value 
as L — > +oo that is significantly larger than the diffraction limit X oc /(29), where A oc is the carrier 



wavelength in the ocean section [0,d] (see Section 3.9.41. 

Let us note that in the case of negligible radiation losses we have the same diffusive continuous 



model as in Theorem 3.7 but with two reflecting boundary conditions at u = and u = 1, (see Section 
4.4 in Appendix to see the full details). However, in this case the solution of the diffusion equation is 
trivial because of the energy conservation. As a result the main difference between the cases with or 
without radiation losses is the low pass filter effect produced by the absorbing boundary condition at 
u = 1 (the bottom of the waveguide), which breaks the energy conservation. 

Let us also remark that the exponential decay of the refocused wave amplitude can also be obtained 



directly from the diffusive continuous model [15, Theorem 6.3] described in Theorem 3.7 Finally, we 
also have the same result for otyi = 1, that is the carrier wavelength and the magnitude of the 
time-reversal mirror are of same order. 

Theorem 3.8 For olm = 1, in the continuum limit N(uiq) 3> 1, we have 

lim H 1 Xg (u Q ,.,L)-H 1 Xo ( C j ,.,L)=0 

in E l , which is the topological dual of E equipped with the weak topology, and where 

di 



lim H 

UJQ— >- + 00 



X, 



x 



6- 



H(x,L). 



Here, H(x, L) is defined in Proposition 3. 1 



Let us remark that in the case of a random waveguide, the order of magnitude ctu of the time- 
reversal mirror plays no role in the transverse profile compared to the homogeneous case. 



3.9.4 Numerical Illustrations 

In this section we illustrate the spatial focusing of the refocused wave around the source location. First, 
we represent the evolution of 7i(L, u), in presence of radiation losses, with respect to L. Here, 7i(L, u) 
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Figure 9: Normalized transverse profile. In (a) and (b) the dashed curves are the transverse profiles in 
the case where the radiation losses are negligible, and the solid curves represent the transverse profile 
H (x, L). In (a) we represent H(x, L) with L = 75, and in (6) we represent H{x 1 L) with L = 250. 



is the mean mode power for the [jV(o;o)w]th propagating mode in the continuum limit N(u>q) 1, 
which is the solution of the partial differential equation in Theorem |3.7 



Second, we represent the transverse profile H(x,L) defined by (34 1 of the rcfocuscd wave, and 
finally we illustrate the resolution of the refocused wave as the propagation distance L becomes large. 

In this section, we consider the following values of the parameters. For the sake of simplicity, we 
take ao = 1, and the correlation length of the random inhomogencitics in the longitudinal direction 
is 1 (a = 1). Moreover, we take ri\ = 2 for index of refraction in the ocean section [0,d], and depth 
d = 20. 



We saw in Theorem 3.7 and Theorem 3.8 that 7i(L, u), in the presence of radiation losses, plays an 
important role in the transverse profile of the refocused wave. In Figure [5J we illustrate the influence 
of the radiation losses on Ti(L,u) as the propagation distance L increase. As we can see in Figure [9] 
and Figure [lOj the radiation losses degrade the quality of the refocusing. Moreover, for L ^> 1, one 
can see a threshold of the quality of the resolution since 



»1 rl 

,AiL 



H Xo (x,L) ~ e 1 / 4>oo,i{v)dv I (j>oo,i(u) cos(2Trxu)du, 
L>1 Jo Jo 

where Xi < is a simple eigenvalue, and the largest one, with corresponding eigenvector ^ j of the 
diffusion operator introduced in Theorem |3.7| (see [TrU Lemma 2.2]). 
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Figure 10: Representation of the evolution of the resolution with respect to the propagation distance 
L. 



Conclusion 

In this paper we have analyzed the time reversal of waves of a broadband pulse in an underwater 
acoustic channel with random perturbations. In this context, using the continuous diffusive models 
developed in |15j . describing the mode-power coupling between the propagating and radiating modes, 
we describe in a simple way the refocused transverse profile (Theorem 3.7 and Theorem 3.8 1 to lighten 



the negative effects of the radiative losses on the time-reversal refocusing property. We have seen that 
radiation losses degrade the quality of the refocused transverse profile as the propagation distance 
increases. First, the amplitude of the refocused wave decays exponentially with the propagation 



distance (Section 3.6 1. Second, using a low pass filter representation of the refocused transverse 
profile (Theorem 3.7), we have shown that the width of the main focal spot increases and converges 
to an asymptotic value as the size of the random section increase, which is significantly larger than 
the diffraction limit X oc /(29) obtained in Propositio n |3.1| (whe re A oc is the carrier wavelength in the 
ocean section with index of refraction m, and 8 = y/1 — l/nf). 



4 Appendix 



4.1 Proof of Proposition 3.1| 

Let us first note that according to Lemma [l.l| and ([7]), we have 

sup \A? - - 1 = 0(N<*-') 

j£{l,...,N-N a } a 



(35) 



where Aj is defined by (|5j. Then, we split the transverse profile in two part in order to use Lemma 
[Til and f35), 



J \l-a M 



d 9 + di 



x 



~e~ x ) = ~2e 



AU 



(f>j(u) ,X + XocX/9)-^ 



1 — cos 



N-[N a ] N 

E + E 

j = l j=N-[N«]+l 

2d M + Kc I {d2-di) 

d 



sine a 



'3 d 

so that the second sum on the right of the previous equality is of order 0(N a ~ 1 ). Moreover, 



N-[N a ] 



E 0j( u} o,xo)<l)j(uo,xo + \ oc x/9) 

3=1 

Ajd 
2 



2d 



±&^y nc ^.&k±M^ | < ^ HN) 
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Now, for the first sum of the previous equality we have 

A 2 



Od 



X — COS 



2x$ + X oc x/9 



and 



/ 2x + X oc x/9\ f . 

i [aj J = cos ^{a 3 - jir) 

- sin ((o-j - j%) 



2xq + X oc x/9 



IS M7T 



2xq + X oc x/9 



2x + X oc x/ 



sin l j7r 



_2x + X oc x/9 



Then, using the Abel transform and Lemma we get 

N-[N a ] 



Xoc\ cosier 



3=1 



< #jV5 



Moreover, using (35 1, we also have 

N-[N" 



X nr ,d 



N-[N° 



3 = 1 



3 = 1 



with 



1™ ^tt; ^ cos f^TT 71 "^ = / cos(2u7ri)du = sinc(27ra) 
o->+oo 2t)d z — ' \ N J In 



3 = 1 



That concludes the proof of Proposition |3.1 



4.2 Proof of Theorem l3?fl 

The refocused transversed profile is given by 



d XI 



d 2 + d x 



-H 



N 



<" + -7rx,L] = ^2 Tj(u} ,L)<j) l (u}o,xo) 



A 2 d 



(f>i(uj ,xo + X oc x/9)—^ 



1 — cos 



29 

3,1=1 

2d M + X^'(d 2 -d 1 ) 
d 



sinc^cr 



X^(d 2 + d 1 ) 
d 



Using the following probabilistic representation 



" Jo A y« (w ) 1 (5- s jv =iv)u(y«=jv-i) ds 1 



Y N = l 



where (Y t N ^ t>Q is a jump Markov process, with state space {1, . . . , N}, intensity matrix T c (ojo), and 
invariant measure /ijv, the uniform distribution over {1, . . . , N}, we have 
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^ T-(uJo,L)(j>i(u!Q,xo)(f>i(uJo,x + X oc x/8) 



AU 



cos 



2d M + X^{d 2 -d 1 ) 



sine a 



Kc M (d 2 + d 



j 



N 1 1 

< Kc" M N [ ^fZJ) ¥ ^ ( Y l = j) + if» N (Y{? = 1) 



3=2 

< KX\- aM hx(N). 
Consequently, the transverse profile of the refocused wave is given by 



~29 



N AU 



^2 -k~Tj (wo, L)(j>i(uj , x )cj)i(ujQ, x + X oc x/9). 



3,1=1 
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Let 77 > such that 77 1 . Thanks to (351, we have 

oc ^ — ^ , 
-20- 2^ ~2~Tj (uo,L)(j)i(uJo,xo)(j)i(u)Q,xo + X oc x/0) 

3,1=1 



[jv(i-n)] 



A r „ ^"'A 2 d 

"26 ^ 2 

3,1=1 

+ o( v ) 

, [JV(i-n)] 

3,1=1 

+ 0( V ). 



Let /♦'(t;) = 1 [0 , i- v ](v), we have 
A 



w AU 



E ^ (^0, 0i (^0, a^o )0i (^0,^0 + X oc x/9) 

3,1=1 



Xoc 

"20 



[Af(l-r,)] 

^2 Tjv (u ,L)(/)i(ujq,x ) <j>i (uj , x + X oc x/9) 
1=1 
0{rj). 



Now, we are able to use the continuous diffusion approximation given in |15l Theorem 6.2]. In what 
follows 7/i is the solution of the diffusion equation given in Theorem 3.7 but with initial conditions 
f v . As a result, we have 

j [N(l-v)] 

- \T},(u> ,L)-T f n{L,l/N) 



1=1 



[N(l-V)]-1 r (l + l)/N 



< 



E 

1=1 



Tj" ul (« 0l L) - T P (L, [Nu}/N) \du 



l/N 



— / \Tfi,(oj ,L,u) —Tf*(L,u)\du 
Jo 



T f v(L,u)- T f v (L, [Nu]/N)\du, 
where the terms on the right side of the last inequality converge to as u>o — > +00. Then 

Xn 



N AU 



■)0 E —^-7~j(uo,L)(f)i(uJo,x )<f>i(uj ,xo + \ oc x/6) 
3,1=1 



V 

' 29 



[iv(i-n)] 



Tfv(L, l/N)<f>i(u ,xo)(f)i(uJo,x + X oc x/9) 



<D{rj). 



Moreover, we have 



A A 2 

(f>j(u ,x )<j)j(u) Q ,XQ+ -j^x) = -± 



2x$ + X oc x/9 



and 



cos (a 



2x + X oc x/9 



) = cos (Oj -3*) 



2xq + X oc x/9 



) cos {jw- 



2xq + X oc x/9 



sin ( (o-j - jit 



2x + X oc xi 



sin l j7r 



_2x + X oc x/9 
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Using the Abel transform, Lemma p5"| , and the continuity of v H> Tfm (L, v) on [0, 1], we get 

\N(l-r,)] „ > 

^lim^A^I £ 7), (£, i/^M? cos (uj ° + ^° c 7 ) | = 0. 
Moreover, using (35 1 on more time and the fact that lim Wo sup ? - \ oc \&j — = 0, 



A IN{1 ~ ,1)] / A \ 

55- E W^/A^cos^) 



A [N{1 ~ V)] / Z x 

= ^2^ E r / ,(L,Z/iV)cos(2.- 5 ) 

= (1 — ? ?) / T}v(L,u) cos(2Trux)du. 



Consequently, according to the proof of |15( Theorem 6.2], we have 

||r / ,(L,.)-Ti(L ) .)l|i»([o, l ]) < WP-MUHm), 

and then 

Xoc ^ AH . 

2^ (wo,i)^(wo,x )(/>z(a;o,xo + X oc x/9) 



lim 

cjq— >+oo 



26> ^ 2 



l 



7i(L, v) cos(2irux)du 



< Krj. 



This concludes the proof of Theorem |3.7| B 

4.3 Proof of Theorem Q 

Let M > 1 and f M = Y,f=i Uj&j e £ M- Moreover, let 



M 



Using (35 1 and 



we have 



Finally, by considering 



sup |<7j -jir\ = 0[ i 
je{i,...,M} 



JV 



sup |/»-/»| = o(i 



we have 

(M\ I TjX (, , r\ £V1 /, , r\ fM fM 



( ff x„( w o, •>£) ~ Hl o {w Ql .,L), f ) L2(0 d) -(i^ (w , .,L) - iJ* ( w o, -,L),f M - / M ) i 2 (0id ) 

+ (Hl {u , .,L) - HI Q {U , ,,L), f ) L 2( jy 



with 



|«(w , ,L) - iT^o, .,L), f M -f M ) LH04) \ 



r N i 

AT ' 



29 



and 



I ( H x (uo,; L)~Hl (ui ,.,L), f M ) L2 (od) \ 



M 

<^E [V^o, L) + p^jF^ € {[Nr,] + 1, . . . , N}\Y N = 



for ?/ > 0, and f v (v) = l[o 1?7 ]- Therefore, it suffices to study J2f=i 7fn ( w 0j £)■ To do this, let g'' be a 
smooth function with compact support included in [0,2^] and such that < f v < g v < f 2r '. Using 
the second part of [T5j Theorem 6.2], 



M 



lim V 7l (w , L) = MT g n (L, 0) < MP (x(i) e [0, 2tj]) . 



1=1 



Here, we have used the probabilistic representation of Tg V (wo,L) introduced in the proof of [T51 
Theorem 6.2], where Po is the unique solution of a martingale problem starting from 0. However, the 
probabilistic representation can be chosen such that the associated diffusion process has aabsolutely 
continuous transition probabilities with respect to the Lebesgue measure Therefore, 



and the rest of the proof is the same as the one of Theorem |3.7[ ■ 



0, 



4.4 Proof of Proposition 3.5 



Following the proof of Theorem 3.7 and using |15[ Theorem 6.4], we get 

"i 



lim 



A: 



H 



0,a M ( 



X, 



uj a ,x 



-x,L) = 



7i(L, u) cos(2ttux)cIu, 



where 7~i(z,v) is a solution of 



d d 
dz T ^ u) = du 



with boundary conditions 

d_ 

du 



Ti(z,0) = 0, 



du 



Ti(z,l)=0, and 7I(0,«)=1. 



However, thanks to the conservation energy caused by the two reflecting boundary conditions, and 
because the initial condition is constant equal to one, this problem admits only one solution, which 
is Tx(z, u) = 1. 

The transverse profile of the refocused wave is studied using the diffusive continuous model in- 
troduced in [T3], with two reflecting boundary conditions at u — (the top of the waveguide) and 
u = 1 (the bottom of the waveguide). Here, the two reflecting boundary conditions mean that there 
is no radiative loss anymore (see Figure |6|, and then the energy is conserved. This is the reason why 
T~x(z, u) — 1. Consequently, the sine profile obtained in Proposition 3.5 is the best transverse profile 
that we can obtain. ■ 
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